Fixed scale approach to the equation of state on the lattice 
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We propose a fixed scale approach to calculate the equation of state (EOS) in lattice QCD. In this 
approach, the temperature T is varied by N, at fixed lattice spacings. This enables us to reduce 
T — simulations which are required to provide basic data in finite temperature studies but are 
quite expensive in the conventional fixed-A^ f approach. Since the conventional integral method 
to obtain the pressure is inapplicable at fixed scale, we introduce a new method, "T-integration 
method", to calculate pressure non-perturbatively. We test the fixed scale approach armed with 
the T-integral method in quenched QCD on isotropic and anisotropic lattices. Our method is 
found to be powerful to obtain reliable results for the equation of state, especially at intermediate 
and low temperatures. Reduction of the computational cost of T — simulations is indispensable 
to study EOS in QCD with dynamical quarks. The status of our study in Nf — 2 + 1 QCD with 
improved Wilson quarks is also reported. 



1. The method 

The equation of state (EOS) of hot QCD is an essential input to understand thermal properties 
of the quark matter in early Universe and in relativistic heavy ion collisions. The lattice QCD is 
the only systematic method to calculate EOS for wide range of T across the phase transition. The 
temperature on the lattice is given by the inverse of the temporal lattice extent: T = (JV f a) _I with 
N, the lattice size in the euclidian temporal direction and a the lattice spacing. Conventionally, 
finite temperature lattice simulations are performed in the fixed-Af, scheme, in which T is varied 
by changing a (or equivalently the lattice gauge coupling [} = 6/g 2 ) at a fixed N,. Because 
this method requires a large amount of computational resources, most studies have been done 
with staggered-type lattice quarks which are less demanding among lattice quarks but whose 
fundamental properties such as locality and universality are not well established. Therefore, to 
evaluate the effects of lattice artifacts, it is indispensable to carry out studies using theoretically 
sound lattice quarks, such as the Wilson-type quarks. For this purpose, we developed the fixed 
scale approach with the T-integral method which can reduce numerical costs 0JJ] . 

To study EOS, it is convenient to first calculate the trace anomaly e - 3p where e is the energy 
density and p is the pressure. For the case of quenched QCD, 
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where a-£ t is the lattice beta-function whose non-perturbative values are relatively easy to eval- 
uate. The subtraction of the zero-temperature counterpart is required for renormalization. To 
obtain a non-perturbative estimate of p, the integral method is usually adopted: Using the ther- 
modynamic relation p = (T/V) InZ valid in the large volume limit, 




with p(J3o) ~ 0. In this conventional method, major part of the computational cost is devoted to 
zero temperature simulations to set the lattice scale, to carry out zero-temperature subtractions, 
and to determine the non-perturbative beta functions, at each simulation point. In the full QCD, 
lines of constant physics (LCP's) should be determined too at T — 0. 

In this report, we push an alternative approach where temperature is varied by N, with other 
parameters fixed, i.e. at a fixed point in the coupling parameter space. While this enables us to 
commonly use one T — simulation for all temperatures, the conventional integral method in the 
coupling parameter space is inapplicable. Therefore, we propose a new method, "the T-integral 
method", in yj] to calculate the pressure non-perturbatively: Using a thermodynamic relation 
valid at vanishing chemical potential, we obtain 




with p(T ) * 0. 

Since the coupling parameters are common to all temperatures, T = subtractions can be 
done by a common zero temperature simulation, the condition to follow the LCP is obviously 
satisfied, and the lattice scale as well as beta functions are required only at the simulation point. 
Then, the computational cost needed for T — simulations is reduced largely. We may even 
borrow results of existing high precision spectrum studies at T = which are public e.g. on the 
International Lattice Data Grid (ILDG) (3J. Because the lattice spacings in spectrum studies are 
usually smaller than those used in conventional fixed-Af, studies around the critical temperature 
T c , for thermodynamic quantities around T c , we can largely reduce the lattice artifacts due to 
large a and/or small N, over the conventional approach. On the other hand, as T increases, N t 
in our approach becomes small and hence the lattice artifact increases. Therefore, our approach 
is not suitable for studying the high T limit. Note that the merits and demerits of our method 
are complement to the conventional method. Our merits around T c may be a good news for 
phenomenological applications, since temperatures relevant at RHIC and LHC will be at most 
up to a few times T c . 



2. Test in quenched QCD 

We test the method in quenched QCD using the standard plaquette gauge action Q}]. Simula- 
tion parameters are summarized in Table [T] The il, i2 and i3 lattices are isotropic, while the a3 
lattice is anisotropic with £ = a s /a t = 4. The temperature ranges cover T ~ 200-700 MeV. 

The trace anomaly obtained on isotropic lattices are summarized in the left panel of FigQ] 
The shaded line represents the result of the conventional fixed-N, method obtained on a quite 
large lattice of N, = 8 and the spatial lattice size N s = 32 (about 2.7 fm around T c ~ 290 MeV) 
[2]. Comparing the results on il and i3 lattices, we find that the lattice cutoff effects are quite 
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Table 1: Simulation parameters on isotropic and anisotropic lattices 1 1]. The beta function is taken from Q|. Corre- 
sponding T = simulations are done on N, = 20 if lattices. 




small on these lattices. On the other hand, the i2 and N t — 8 lattices show small deviations near 
T c . These deviations may be explained by the physical finite size effect expected in the critical 
region. Off the critical region, all results agree well with each other. 

Dotted lines in the left panel of Figfflare the natural cubic spline interpolations. We have 
sufficiently many data to smoothly interpolate. Carrying out the numerical integration OJ, we 
obtain the pressure shown in the right panel. Error bars shown in the figures represent the sta- 
tistical errors estimated by the jackknife method. To estimate the systematic error due to the 
interpolation ansatz, we show the results of interpolations with the trapezoidal rule (labeled by 
"T") too. We find that all the results agree well with each other within the statistical errors. 

In the T-integral method, T is restricted to discrete values due to the discreteness of N t . We 




attempt to directly estimate systematic effects due to discreteness of T by comparing the results 
with those obtained on an anisotropic lattice a2, which has about 4 times finer resolution in T 
than the i2 lattice. In Figf2]left, we compare the trace anomaly on a2 and i2 lattices. We find 
that the data points from the a2 lattice appear on the interpolation line from the i2 lattice, while, 
due to the cruder resolution in T, the natural cubic spline interpolation for the i2 lattice slightly 
overshoots the data on the a2 lattice around the peak. We note that the height of the peak on 
the a2 lattice is consistent to those of the fine i3 and N t - % lattices. Carrying out the numerical 
integration, we obtain EOS shown in Figf2]right. We find that the results are well consistent with 
each other, suggesting that the systematic errors due to the discreteness of T is at most about the 
statistical errors. See Ref.[ 1] for more discussions. 

3. Outlook towards the EOS in QCD with dynamical u, d, s quarks 

We have presented a new way to calculate EOS on the lattice. Our method is complemen- 
tary to the conventional method, and will have an advantage around the transition temperature. 
We note that our method can be easily applied to finite densities using the conventional Taylor 
expansion method. 

We are currently investigating EOS in 2 + 1 flavor QCD with non-perturbatively improved 
Wilson quarks, using the configurations by the CP-PACS/JLQCD Collaboration [4], which are 
public on the ILDG [3]. Among their simulation points, we have chosen the finest lattice (a = 
0.07 fm) with the lightest u and d quarks (m n lm p — 0.63). Using the same coupling parameters, 
we are generating finite temperature configurations on 32 3 x N, lattices with N, = 4, 6, • • •, 16, 
where the pseudo-critical temperature is expected to be around 14. Preliminary results for heavy 
quark free energies on these finite temperature configurations are presented by Maezawa in these 
proceedings 05J, in which another good feature of the fixed scale approach is pointed out: We 
can study the thermal effects on the heavy quark free energies without artificially adjusting the 
constant term of the potential. We hope that we present new results of EOS at the next QM 
symposium. Our final goal is to study EOS in 2 + 1 flavor QCD just at the physical quark 
mass point. We are planning to extend the study using zero-temperature configurations being 
generated by the PACS-CS Collaboration just at the physical point jfj. 
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